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^ ; ABSTRACT 

, We analyze the 2D correlation function of the SDSS-III Baryon Oscillation Spectro- 

ff^ ■ scopic Survey (BOSS) CMASS sample of massive galaxies of the ninth data release to 

measure cosmic expansion H and the angular diameter distance Da at a mean red- 
shift of (z) = 0.57. We apply, for the first time, a new correlation function technique 
called clustering wedges S,Afiis). Using a physically motivated model, the anisotropic 
I baryonic acoustic feature in the galaxy sample is detected at a significance level of 

■ 4.7cr compared to a featureless model. The baryonic acoustic feature is used to ob- 

tain model independent constraints cz/H/rg = 12.28 ± 0.82 (6.7% accuracy) and 
Da/ts = 9.05 ± 0.27 (3.0%) with a correlation coefficient of —0.5, where is the 
sound horizon scale at the end of the baryonic drag era. We conduct thorough tests 
on the data and 600 simulated realizations, finding robustness of the results regard- 
less of the details of the analysis method. Combining with Ts constraints from the 
Cosmic Microwave Background we obtain H{0.57) ~ 90.8 ± 6.2 kms~^Mpc~^ and 
Da (0.57) = 1386 ± 45 Mpc. We use simulations to forecast results of the final BOSS 
CMASS data set. We apply the reconstruction technique on the simulations demon- 
strating that the sharpening of the anisotropic baryonic acoustic feature should im- 
prove the detection as well as tighten constraints of H and Da by ^ 30% on average. 

Key words: cosmological parameters, large scale structure of the universe, distance 
scale 
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One of the most exciting recent observations is the acceler- 
ation of the expansion of the Universe since redshift z = 1 
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(iRiess et al ] |l998l : IPerlmutter et al.|[l999l l . The origin of this 
phenomenon is thought to be an energy component with 
negative pressure (e.g, the so-called dark energy or a cos- 
mological constant), or o therwise a break-down of General 
Relativity ()Einsteinlll916h on cosmic scales. For an in-depth 
summary of th e observed acceleration and its possible inter- 
pretations, see IWeinberg etHI \2Ql± . 

One method of measuring geometry from a 3D map of 
cosmological objects is thr ough a geometric technique called 
the Alcock-Paczynski test (|Alcock &l Paczvnskilll97^ ) along 
with a standard ruler k nown as the baryonic acoustic fea- 
ture. [Akoc^^&l^aczxnski (1979) demonstrated that by as- 
suming an incorrect cosmology when converting observed 
redshifts Zobs to comoving distances a spherical cosmo- 
logical body will appear deformed due to geometrical dis- 
tortions. In the context of galaxy maps, this would cause a 
coherent distortion of the apparent radial positions and an 
anisotropic signature in clustering probes. Reproducing an 
isotropic clustering signal would result in obtaining the true 
cosmology. The observable of this process is HDa, where 
H is the expansion factor and Da is the angular diameter 
distance. To break this degeneracy, and hence improve con- 
straining power, one can apply the technique on a standard 
ruler, such as the baryonic acoustic feature. 

Early Universe plasma-photon waves propagated at 
close to the speed of sound from over dense regions, and 
came to a near halt at the era of decoupling of photons 
from baryons at z, ~ 1100 at a characteristic comoving dis- 
tance of Vs ~ 150 Mpc from the originating over-density. 
This process left a distinctive signature in CM B anisotropies 
and in the large - scale structure of galaxies (|Peebles fc Yul 
Il970l ). IHu et al.l (jl997l ) review how the CMB anisotropies 
can be used to constrain fund amental cosmological parame- 
ters. iBassett fc Hlozekl (|2010l ) review the baryonic acoustic 
signature in the clustering of matter and its usage as a stan- 
dard ruler. 

Following first baryonic acoust ic feature measur e ments 
in t he clustering of galaxies by lEisenstein et al.l (|2005h 
and I Cole et al] (|2005l ). two recent surveys, the W iggleZ 
Dark Energy Sky Survey (jDrinkwater et al] l2010l) . and 



the Sloan Digital Sky Survey (SDSS-III: lYork et al.llioool : 
lEisenstein et al.l 201ll) Barvonic O scillation Spectroscopic 
Survey (BOSS; iDawson et al.|[201^ ). have rep orted detec- 
tions o f t he baryonic acoustic feature at z > 0. 5 ( Blake et al.l 
l2011bl ldl: lAnderson et aLll2012l;[sinchez et al.llioi^ ). as well 
as the 6dF Galaxy S urvev (| Jones et"al] 20091 ) at z < 0.2 
Beutler et al] l201ll ). iBusca et al.l (|20ia ) and ISlosar et ahl 
20131 ) also detect the baryonic acoustic feature, for the first 



time, in the Lyman-alpha forest of BOSS quasars between 
2 < 2 < 3. 

The focus of of most of these studies has been on the 
angle-averaged signal which constrains (D\/ H)'S /va, where 
rg is the sound horizon at the end of the baryon drag era 
(see !j2]). This degeneracy originates because for every line- 
of-sight clustering mode (which constrains Hr^), there are 
two transverse modes that constrain Da/tb- 

The subject of this study is breaking the D\/H de- 
generacy by using the Alcock-Paczynski effect through 
anisotropic clu s tering . This approach was first suggested by 
IHu fc Haimaiil (|2003|) bv usin g the two-dimensional power 
spectrum P(k). 1wagner et al.l (|2008l ) used mock catalogues 
at z = 1 and 3 to demonstrate the usefulness of the tech- 



nique. IShoii et al.l (|2009l ) argued that H and Da information 
is encoded in the full 2D shape, and presented a generic al- 
gorithm that takes into account dynamic distortions on all 
scales, assuming all non-linear effects are understood. First 
attempts to apply these technique s on 2D P(k) and ^(s ) 
clustering planes were performed bv Okumura et al] (|2008h . 
IChuan g & Wang (2011), and Bla ke et al.l l|201ld ). 

Padmanabhan fc Whitel (|200^ ) suggested decomposing 



the 2D correlation function into Legendre moments. They 
argue that the monopole (^o angle averaged signal) and 
the quadrupole components (^2, see Equatio n 1101) contain 
most of the relevant constrain ing information. iTaruva et al.l 
(|201ll ) and lKazin eFal] (|2012h show that the hexadecapole 
term ^4 contains extra constraining power, which could be 
harnessed in the future with higher S /N than that currently 
available. 

The advantage of analyzing ID projections over the 2D 
plane is the relative simplicity of building a stable covariance 
matrix. 

The first analyses of the anisotropic baryonic acoustic 
feature using ^0 and ^2 have been perform ed on the SDSS- 
II luminous red galaxy sample {z ~ 0.35; IXu et aL I l2012al : 
IChuang fc Wand l2012allbl) and the p R9-CMASS sample 
tested here {z ~ 0.57: iReid et al.ll2012l ). 

We analyz e, for the fir s t tim e, an alternative ID ba- 
sis sugg ested by Kazin et al.l (|2012| ) , called clustering wedges 
f"A,i(s). rGaztanaga et al. ( 20091 ) focused on a narrow cluster- 
ing cylinder ^(sh , s j_ < 5 /i~^Mpc). In a subsequent analysis, 
iKazin et al.l (|2010l ) proposed using wider clustering wedges 
f Au(s ) to improve S/N of the measurements, [kazin e t al] 
( 2012h analyzed the constraining power of H and Da of 
^Aji=o.5(s) on mock catalogues. They concluded that these 
statistics should be comparable in performance to the mul- 
tipoles (^0, C2) and provide a useful tool to test for system- 
atics. The current study is the first analysis to perform such 
a thorough comparison on both data and mock galaxy cat- 
alogs. 

Our analysis differs from the previous ones in a few 
other aspects. First, we compare results both before and af- 
ter reconstruction. Reconstruction is a technique which cor- 
rects for the damping of the baryonic acoustic feature due 
to the large-scale coherent motions of galaxies. The bary- 
onic acoustic feature is sharpened by calculating the dis- 
placement fi^ld_aiid_sluftirig^ to their near-original 
positions (|Eisenstein et al.' 2007ah. Second, we follow a sim- 
ilar approach as in Xu et al. ( 2012al ). by focusing on cz/ H /r^ 
and Da/ts and margin alizing over shape information. One 
notable difference from IXu et al.l (l2012al ). however, is that 
they apply a linear approximation of the Alcock-Paczynski 
test, where here we use the full non- linear equations. We 
compare both methods in Appendix [EJ and show that the 
linear approach under-estimates the uncertainties of the ob- 
tained constraints. Finally, we compare between two inde- 
pendent theoretical ^ templates. 

This study is part of a series of papers analyzing the 
anisotropic clustering signal of the DR9 CMASS galaxy sam- 
ple, containing 264, 283 massive galaxies between 0.43 < z < 
0.7. Here we measure H and Da in a mo del independent 
fashion through = 0.5 clustering wedges. [Anderson et al.l 
(|2013h uses "consensus" values of clustering wedges and mul- 
tipoles to infer cosmological implications. Both of these stud- 
ies focus on the information contained within the anisotropic 
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baryonic acoustic feature. Two furthe r studies analyze the 
information from the full shape of ^ f slilSanchez et al] (|2013l ) 



use the 5ah=o.5, and lChuang et a~ (|2013h focus on the : 
tipoles ^0,2. 

This study is constructed as follows: in iJ2]we explain in 
detail the geometric information encoded in redshift maps. 
In !j3]we define the clustering wedges and in ^we present 
the data and mock catalogs. In !jS]we describe the method 
used in our analysis; ij6]describes our results. We discuss the 
results in iJ7]and summarize in [JS] 

To avoid semantic confusion, we briefly explain here the 
terminology of the different spaces mentioned throughout 
the text. First, all analyses are based on two-point corre- 
lation functions, which we refer to as configuration-space, 
as opposed to the Fourier domain called fc— space. Second, 
when referring to a space affected by redshift distortions, we 
call it redshift space, and when there are none we refer to it 
as real space. 

All the fiducial values calculated her e are based on 
using the WMAP7 flat ACDM cosmology (|Komatsu et al.l 
l201lf ). To calculate comoving distances we assume the mat- 
ter density ilu = 0.274. Assuming h = 0.7 this yields 
fiducial values: = 93.57 kms-^Mpc-\ = 1359.6 



Mpc at 



0.57. Throughout, we also use derived unit- 



less relationships (cz/H/rsY = 11.94, (DA/rs) = 8.88, 

: 1020). For these we as- 
0.0224, radiation density 



where rl = 153. 1 Mpc (at 
sume the baryon density Qth^ 



lO^Qrh^ = 4.17 and photon density of lO^n^h'-^ = 2.47. 



2 COSMIC GEOMETRY FROM GALAXY 
MAPS 

Although galaxy distributions in real-space are assumed to 
be statistically isotropic, measured clustering signals from 
galax;ies from redshift maps are anisotropic. This is a re- 
sult of two physical effects that are at play when converting 
observed redshifts Zohs to comoving distances x'- 



XiZob 



dz 



(1) 



The first, which we refer to as redshift-distortions, stems 
from the fact that ^obs is a degenerate combination of the 
cosmological flow and the radial component of the peculiar 
velocity. This results in anisotrop ic clustering components 
due to large-scale coherent flows (|Kaisei|[l987l ). and veloc- 
ity dispersion effects within galaxy clusters. For a detailed 
introd uction on dynamical redshift-distortions see lHamiltonI 
(|l998h . 

On large scales, these effects can be used t o test for de- 
viat ions from General Relativity ( Kaiser! Il987l: iLindeil 20081: 
see IGuzzo ct al.' '2007*; 'Samushia ct al."20I2!: iBlake et all 
2011a : Samushia ct al..,2013 : Bcutlcr ct al..,201 j for the most 
recent measurements). The observable in this test is fas, 
where b is the linear tracer to matter density bias, / = 
dDi /d In a is the rate of change of growth of structure, Di is 
the linear growth of structure, and erg is the linear r.m.s of 
density fluctuations averaged in spheres of radii 8/i~^Mpc. 
This study focuses on a second more subtle effect which in- 
volves geometric distortions. 

Comoving separations between two nearby points in 
space depend both on z and the observer angle between 



them 0. Assuming the plane-parallel approximation be- 
tween galaxy pairs, radial separations are defined as S|| = 
cAz/H{z), where c is the speed of light, and transverse dis- 
tances s± = 0(1 + z)Da, where the proper (physical) angu- 
lar diameter distance Da is defined as: 



Da 



1 



1 



l + z Ho 



X 



c/Ho 



(2) 



where Hq = H(0) and SIk = 1 — J^x is the represen- 
tation of the curvature, and Qx are the energy densities of 
compoenents X (matter, radiation, etc.) Q Hence, assuming 
an incorrect cosmology in Equation Q would cause a spher- 
ical body (meaning S|| = s±) to be deformed. For example, 
a lower H{z) than the true one would cause an elongation 
along the line-of-sight due to an increased S|| , where a lower 
Da{z) than the true value would cause a transverse squash- 
ing, because of a decrease of s± . Therefore, by fixing the ob- 
servables Q and Az, retrieving a spherical shape constrains 
the HDa combination. 

Various techniques have been suggested to mea- 
sure HDa through this Alcock-Pac z ynski tes t (AP 
henceforth; lAlcock fc Paczvnskil 1 19791 : IPhillippd 1 1994 
iLavaux fc Wandeltl |2012| ). Here we focus on clustering of 
galaxies, where line-of-sight clustering modes depend on S| 
(1/H) and transverse modes on s± (Da), and hence the 
anisotropics due the AP effect. 

It has been pointed out that the anisotropies from this 
geometr ic effect are de g enera te with those from redshift- 
effects |Ba llingcr ot al. ' 119961). Various studies, such as 
iBlake et al. (2011c) and Rcid ct al . (20l2), show the degen- 
eracy between HDa and fag. In this study we marginalize 
over the redshift distortion information and focus on the 
geometric distortions. 

In practice, when converting redshifts to comoving 
distances, the Ho factors out trivially and thus we ex- 
press comoving distance in units of h~^Mpc, where h = 
//o/(100kms~^Mpc~^). The rest of the parameters in H{z) 
(Qx and their equation of states wx) have more important, 
and potentially measurable, effects. 

One way of overcoming these effects is to recalculate x 
and then the clustering statistics for every set of parameters 
when determining cosmological constraints. However, that 
approach is currently not practical; instead, we vary a fixed 
clustering template, as described below. 

Although the baryonic acoustic feature comoving scale 
is fixed, the apparent position measured in the correla- 
tio n function depen d s on H ts and Da/ts- As demonstrated 
in lEisenstein et al.l (120051 ). the baryonic acoustic feature 
in the angle average signal is sensitive, to first order, to 
{D\/hY'^ /ts. iPadmanabhan fc Whit j (|2008l ) show that 
analysis of the anistropic signal adds HDa information, and 
hence breaks the degeneracy. To break the degeneracy with 
Ts one needs to add additional information from the CMB 
anisotropies. 

When relating measured from the CMB to that in the 
large-scale structure, one must take into account that these 
two definitions correspo nd to slightly dif ferent sound horizon 
radii (see Equation 1 in lBlake fc Glazei3rook.200 j ). Because 
the baryons have momentum at decoupling z«, the baryonic 



Note that this is generic because isin(ix) = — sinh(a;). 
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acoustic signature in the distribution of matter is related to 
i"s(2d) > rs{Zf), where is the epoch when the baryonic 
drag effectively ended jEisenstein fc Hul Il998l ) . The bary- 
onic acoustic signature in the CMB anisotrop ies corresponds 
to Zf For current r^i^Zf) measurements see iHinshaw et all 
(|2012[). and for r^iz d ) pred ictions from the CMB, see Table 
3 of lKomatsu et al.l lj2009l l. 

Conservation of the observer angle means that true 
separations transverse to the line-of-sight component s\ will 
be related to an apparent "fiducial" component s'l by|j 



■ ax, 



where 



(3) 



(4) 



where the "f ' subscript indicates the fiducial cosmology when 
calculating x{z), and "t" indicates the true cosmology. 

Similarly, the true line-of-sight separation component is 
related to the fiducial by: 



4, 



with 



^11' 



(5) 



(6) 



The sound horizon r-s(zd) terms appear due to the degener- 
acy with Da and H, when applied to the baryonic acoustic 
feature as a standard ruler. Here we quote the rescaling in 
the position of the peak of the ^. The purely geometrical 
effect of changing the cosmology does not depend on rs(2d). 

In Appendix|X]we explain how we apply the AP test in 
practice through the mapping of ^ between these coordinates 
systems. 

We also make use of an alternative representation of 
Q|i and a\\ through the isotropic dilation parameter a 



lEisenstein et al.] 2005h and the aniso tropic warping param- 
eter e (|Padmanabhan fc Whit"3l2008l ): 



Da 



H 



DaH 



2/3 1/3 



1/3 



(7) 



(8) 



3 CLUSTERING WEDGES 

Assuming azimuthal statistical symmetry around the line- 
of-sigh10 the 3D correlation function ^(s) can be projected 
into 2D polar coordinates: the comoving separation s and 
the cosine of the angle from the line-of-sight fi, where the 
line-of-sight direction is fi — 1. 

The 2D plane of ^{fi,s) can then be projected to clus- 
tering wedges Afi as: 

CAM(s) = ^y (9) 



^ Here we assume the plane-parallel approximation for each pair. 
^ The assumption of azimuthal statistical symmetry around the 
line-of-sight is true even with geometrical distortions. 



For the purpose of this study, we focus on two clustering 
wedges of — 0.5, which we call line-of-sight (,\\{s) = 
^o.5(Mmin = 0.5, s) and transverse ^±{s) = (fo.5(Atmin = 0, s). 
For consistency we compare all results to the multipole 
statistics defined as: 



= 



21+1 



5(^,s)££(/i)d^, 



(10) 



where Ce{x) are the standard Legendre polynomials. 

The clustering wedges and multipoles are complemen- 
tary b ases of similar information. As shown bv lKazin et al.l 
(|2012h up to order ^ = 4 they are related by: 



1 ^ 



a(s) = eo(s)- 



(11) 



(12) 



A useful relationship is the fact that the average of the 
Ajj. = 0.5 clustering wedges results in ^o- 

In real space, where there are no anisotropies, all ^ > 
components are nulled, and clustering wedges of any 
width correspond to the monopole signal|f| The AP effect 
breaks this symmetry, causing £ > components due to 
geometric distortions. 



4 DATA 

We base our measurements of cz/H/va and Da/ts on the 
large-scale anisotropic correlation function of the BOSS 
DR9-CMASS galaxy sample. Here we give a brief description 
of the sample, and the calculated ^(/i, s). 

4.1 The DR9-CMASS galaxy sample 

We use data from the SDSS- III BOSS survey 

(lEisenstein et al. 1 1201 ll : iDawson et all |2013| ). The galaxy 
targets for BOSS are divided into two samples, LOWZ and 
CM ASS. These are selected on the basis of photometric 
observat ions carried out w i th a d rift-scanning mosaic CCD 
camera jGunnet al.lll998l . l2006h on the Sloan Foundation 
telescope at the Apache Point observatory. Spectra of 
these galaxies are obtained using the double-armed BOSS 
spectrographs (|Smee et al.l |2012| ). Spectroscopic redshifts 
are then measured by means of the minimum - y^ tem plate- 
fitting procedure described in lAihara et al.1 (120111 ) ■ with 
te mplates and method s updated for BOSS data as described 
in lBohon et~aLl (|2012l '). 

Our analysis is based on the CMASS galaxy sam- 
ple of SDSS Data Release 9 (DR9. lAhn et al.ll2012l ). This 
sample was designed to cover the redshift ranges 0.43 < 
z < 0.7 down to a limiting stellar mass, resulting on a 
roughly constant number density of n ~ 3 x 10~'*/i'^Mpc^'^ 
(|Eisenstein et al.1 1201 ll : |d awson et al.ll2013l . Padmanabhan 
et al. in preparation). This sample cont ains mostly central 
galaxies, with a ~ 10% satellite fraction (|White et al.ll201ll : 
iNuza et al.]|2012D and it is dominated by early type galaxies, 
althoug h it contains a signifi cant fraction of massive spirals 
(~ 26%. [Masters et al.ll201lh . 

^ Homogeneity and isotropy are assumed here. 
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lAnderson et all (j2012l ) present a detailed description of 
the construction of a CMASS catalogue for LSS studies. We 
follow the procedure detailed there and refer the reader to 
this article for more details. 

4.2 PTHalo mock catalogues 

Mock catalogues play a major role in the analysis and inter- 
pretation of large-scale structure information, as they offer 
a useful tool to test for systematics and provide the means 
with which to estimate statistical errors. In this analysis 
we use 600 PTHalo mock galaxy realizations to test our 
analysis pipeline and construct a covariance matrix of our 
m easurements. F u ll deta ils of the mock catalogues are given 
m iManera et al] (|2012h . Briefly, the mocks are based on 
dark matter 2LPT simulations (2""^ order Lagrangian Per- 
turbation Theory) , that were populated with mock galax- 
ies w ithin dark matter haloe s . The halo occupation distribu- 
tion jPeacock fc Smithll2000l: IScliak"2000': IScoccimarro et al.l 



l200l| - iBerlind fc Weinberg| |2002: Coorav fc ShethH2002l ') pa- 



rameters are determined by comparing the correlation func- 
tion to that of the data in the scale range of [30, 80] /i^^Mpc. 

To match the selection function of the data, the mock 
data is split by the Northern and Southern CMASS angular 
geometry and galaxies were excluded to match the radial 
profile. 

4.3 The Anisotropic Correlation Function 

To compute the co rrelation function, we use the 
iLandv fc Szalavl () 19931 ) estimator with an angular depen- 
dence: 



DD(fi, s) + RR{fj., s) - 2DR(fi, s) 



(13) 



RR{fi,s) 

We calculate the normalized data-data pair counts in bins 
of evenly separated fi and s, DD{fi,s), and similarly for 
data-random pairs, DR, and random-random, RR, where 
for each pair = 1 is defined as the direction in which a 
vector from the observer bisects s. The /i values of each bin 
are the flat mean value. Our choice of binning is A/i — 1/100 
and As = 4/i"^Mpc. 

To obtain the clustering wedges we use Equation ([9]), 
where for ^±{s) we use the jj, range [0,0.5] and for £,\\{s) 
[0.5,1]. The resulting pre-reconstruction clustering wedges 
and multipoles are presented in top and bottom panels of 
Figure[T] respectively. The line-of-sight wedge C!|(a' > 0.5, s) 
is clearly weaker than the transverse wedge ^x(/i < 0.5, s). 
This large difference in amplitudes on large-scales is due to 
redshift-distortions. 

For comparison in the right panels of Figure [T] we show 
the mock-mean signals, i.e., the mean ^a/j and of 600 
mock catalogs. 

4.4 Reconstructing the baryonic acoustic feature 

lEisenstein et al.l (|2007al ) showed that large-scale coherent 
motions, which cause a damping of the baryonic acoustic 
feature, can be ameliorated by using the gravitational po- 
tential estimated from the large-scale galaxy distribution 
to predict the bulk flows, and undo their non-linear ef- 
fect on the density field. First studies focusing on periodic 



boxes shown that this reconstruction technique sharpens 
the baryonic acoustic feature, and hence improves its usage 
as a standard ruler ( Padmanabhan c t al. 1 12OO9I : iNoh et al.l 
l2009l:ISeo eraIi,2010,:,Mchta ct al.,20liirWe follow the pro- 



cedure in Padmanabhan et al. 

I (|2012l ). which takes into ac- 



count practical issue s as edge effects by applying a Weiner 
filter (|koffman fc Ri bak 1991): IZaroubi et al.lll995l '). We ap- 
ply the reconstruction procedure on both the DR9-CMASS 
data, as well as on the mocks. 

Figure [2] displays the post-reconstruction results for 
^ll_x(s) (top left) and the ^0,2(5) (bottom left). We clearly 
see that the amplitudes of the clustering wedges are aligned 
at the scales of the baryonic acoustic feature and larger. 
This is due to the fact that reconstruction not only corrects 
for large-scale coherent motions, but also corrects, to a cer- 
tain extent, for redshift-distortions, as is seen by the near 
nullifying of the ^2{s). 

For comparison, the right panels Figure [2] show results 
of the post-reconstruction mock-mean signal. We clearly see 
that the ^2(s) reverses from negative at baryonic acoustic 
feature scales from the pre-reconstruction signal to positive. 
This change might be attributed to an over compensation 
of the redshift-distortions. In other words, throughout the 
reconstruction process, we estimate / to shift galaxies in 
the radial direction, with the aim to reduce the Kaiser effect. 
An over-estimation could potentially put field galaxies a bit 
further away from high dense regions, and hence reverse the 
^2(3) signal, yielding a £,±{s) that is slightly weaker than the 
5ll(s). We are not concerned with this issue, because we do 
not expect redshift-distortions to shift the position of the 
anisotropic baryonic acoustic feature. 

In both clustering wedges, in pre- and post- 
reconstruction, there is a clear signature of the baryonic 
acoustic feature. We quantify the significance of the detec- 
tion in i]6.1l 



5 ANALYSIS METHODOLOGY 
5.1 Statistics used 

When computing likelihoods of a model M with a variable 
parameter space $ to fit data D, we calculate the x^'- 



(14) 



where i, j are the bins tested. The likelihood is then assumed 
to be Gaussian -L($) oc exp (— |x^('^*)) ■ 

Throughout this analysis we run Monte Carlo Markov 
Chains (MCMC) nominally for nine or ten parameters as 
described in H5.3I We quote the mode of the posterior as our 
measurement and half the 68% CL region (68CLr hence- 
forth) for the uncertainty, because these are well defined 
regardless of asymmetries in likelihood profiles. 



5.1.1 Covariance matrix 

We construct the covariance matrix dj from the A'mocks ~ 
600 mock catalogues. (For a description of the mocks used 
see H4.2I ) The signals are not independent but have 

significant cross-correlations. To take these correlations into 
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Figure 1. Top Left: the pro-reconstruction DR9-CMASS {z) = 0.57 clustering wedges are displayed multiplied by (s/rs) in the main 
plot, and without in the inset. Bottom Left; the CMASS monopole and quadrupole. The solid lines in the left hand panels are the 
RPT-based best fit models (Renormalized Perturbation Theory; see i)5.2|l . The a-^d degrees of freedom are indicated. Right: The same 
statistics using the mean signal of 600 PTHalo mock realizations. The uncertainty estimates are the square root of the diagonal elements 
of the covariance matrix. The solid lines in the right hand panels are the RPT-based templates. (For clarity the transverse wedge and 
quadrupole are shifted by 1 /i~^Mpc.) The line-of-sight wedge {fi > 0.5 red circles) is clearly weaker than the transverse wedge (fi < 0.5 
blue squares). In both clustering wedges there is a clear signature of the baryonic acoustic feature. 



account, when constructing the dj, we treat the mocks sig- 
nals in an array of the form ^^2b] = [C| | ) 5x] i meaning a ID ar- 
ray with twice the length of the separation range of analysis. 
When analyzing the multipoles we apply a similar conven- 
tion ^[2b] = [CcCa]- We then construct a covariance matrix 
of ^i2a] defined as: 



^ "mocka 



(15) 



Figure [3] shows the correlation matrix Cij / y^d^iCjj 
of C||,x PrS" 3-iid post-reconstrucion. The ^± quartile has 
slightly larger (normalized) ofF-diagonal terms than in the ^| | 
quartile, demonstrated by the less steep gradient. There are 
also non-trivial positive and negative covariance cross-terms 
between the ^± and In the reconstructed we notice a 
sharper gradient, and a shallower negative region, indicating 
less dominance of the off-diagonal terms. This means that 



the reconstruction procedure reduces the covariance between 
the data points. Examining dj of the ^0,2 we find similar 
trends. 

Figure U displays the square root of the diagonal el- 
ements. It is clear that pre-reconstruction the scatter in 
the two clustering wedges is slightly different, where post- 
reconstruction they are similar, and less than that of pre- 
reconstruction. We clearly see that ^0 yields the lowest scat- 
ter of all the ^ statistics, and ^2 the highest, both pre- and 
post-reconstruction. 

To correct for the bias due to the finite number of re- 
alizations used to estimate the covariance matrix and avoid 
underestimation of the parameter constraining region, after 
inverting the matrix to C~,?^^;„,i, we multipl y it by the cor- 
rection factors given in iHartlap et all (|2007l ): 



C"' = d 



-1 

original 



(iVmoc 



A^bins - 2) 



(iVn 



-1) 



(16) 
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Figure 2. Top Left: the post-reconstruction DR9-CMASS (z) = 0.57 clustering wedges are displayed multiplied by (s/rg) in the main 
plot, and without in the inset. Bottom Left: the CMASS monopole and quadrupole. The solid lines in the left hand panels are the 
RPT-based best fit models (Renormalized Perturbation Theory; see i)5.2|l . The a-nci degrees of freedom are indicated. Right: The same 
statistics using the mean signal of 600 PTHalo mock realizations. The uncertainty estimates are the square root of the diagonal elements 
of the covariance matrix. The solid lines in the right hand panels are the RPT-based templates, where we set ^2 = 0, as explained in 
i]5.2.2l (For clarity the transverse wedge and quadrupole arc shifted by 1 /i~^Mpc.) The linc-of-sight wedge {/i > 0.5 red circles) is similar 
to the transverse wedge {/i < 0.5 blue squares). This result shows that reconstruction substantially reduces effects of rcdshift-distortions. 
In both clustering wedges there is a clear signature of the baryonic acoustic feature. 



In our analysis A^mocks ~ 600, A'^bins ~ 76 (2x38) (when 
analyzing region [50,200] h~^Mpc), yielding a factor of 0.87. 



5.2 Non-linear ^ templates 

The modeling is split into two parts: inclusion of redshift- 
distortions and modeling for non-linearities. Here we de- 
scribe the former, and later consider two procedures for 
defining non-linearities. 

Once the non-linear Pnl is defined (see i]5.2.ip , redshift- 
distortions are added such that the non-linear z-space power 
spectrum is: 



PNL{k,fJ.k) 



1 



(l + (fc/avMO')' 



{l + Pd)'PNL, (17) 



where P = f/b. 

Although the velocity dispersion paramet er cry appears 
to be an unresolved subject of investigation (|Taruva et al.l 



I2OO9I ). we find that applying it in the above Lorentzian 
format yields a good agreement with the mock-mean sig- 
nals 5||_x and ^0,2 down to s > 50/i~^Mpc. We find the 
Lorentzian format is preferred over the popular Gaussian. 



The conversion to configuration space is accomplished 
by m eans of Equations (4.8) and (4.17) in iTaruva et al.l 
(|2009l ). As described in Appendix |XJ we apply this calcu- 
lation once to obtain the ^0, ?2 templates, which are stored 
during the MCMC calculations. This approach means that 
we fix the parameters /, /3, crv constant, and allow for their 
effective changes through the ao stat and A{s) shape param- 
eters, as described in >15.3.1I tl6.3.3l The values assumed for 
these parameters are summarized in Table [1] 
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Figure 3. Wo use a suite of 600 PTHalo pre-reconstruction moclc catalogs (Left) and post-reconstruction (Rigiit) to construct tlie 
covariance matrix of tiic clustering wedges, displayed here in correlation matrix form Ci^j j ^fCiJCj2- "^^"^ bottom left quadrant is that 
of the top right quadrant is that of ^x- The other quadrants, which are mirrored, are the cross-correlation between the bins of ^|| 
and 




50 100 150 200 50 100 150 200 



s \hr^ Mpc] s \hr^ Mpc] 

Figure 4. The y/Cu values constructed from the pre- (solid) and post-reconstruction (dashed) mocks. Left plot: results of the clustering 
wedges ^|| (thick red), (thin blue). Right plot: similar for the 6 (thick rod), ^2 (thin blue). Reconstruction substantially reduces the 
covariance in these measurements. 



Table 1. Non linear anisotropic ^ templates 



Template Name 


Equation Base 


Fixed parameter values 






Comment 


^RPT-basod pi-e-rec 


m 


fcNL = 0.19/iMpc~\ Ayic = 


2.44, o-y = 


5.26/i-iMpc 




^RPT-bascd post-rec 


GSli 


fcNL = 0.50/iMpc-\ Ayic = 


2.44, fry = 


Ofe-iMpc 


6 = 


^dewiggled pi-e-rec 


12011 


S|[ = llh.-iMpc, Sx = 6/i" 


"^Mpc, cry 


= Ih-^Mpc 




^dewiggled post-rec 




S|[ = Ex = Sft-^Mpc, (TV = 


= Ih-iMpc 




6 7^ but small 



The RPT-based templates are shown in Figures [T]andj2] 

The dewiggled templates are shown in Figure 1 of lAnderson et al.l feOial V 

After the base equation is calculated. Equation I I17I I includes rodshift-distortions (post-rec assumes /3 = 0) 
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5.2.1 Non-linear P{k) 

We use two anisotropic templates. The primary focus is on a 
physically motivated model based on Renormalized Pertur- 
bation Theory (RPT), which takes into account first order 
corrections of k mode coupling. Throughout this study we 
compare performance of this template to one that describes 
the effect of non-linearities in the baryonic acoustic feature 
through the "dewiggling" procedure. Both templates assume 
an exponential damping of the baryonic acoustic feature due 
to large-scale coherent motions, where in RPT-based this is 
assumed to be isotropic and in dewiggled anisotropic. 
For the RPT-based template we write: 



-PRPT(fc) = PLinoar(fc) CXp ( - ( j j + AmC f Hoop (fc) , 

(18) 



where 



^lloop(fc) - ^ 



dq|F2(k-q,q)|2pLincar(|k-q|)PLin 



ar(q). 

(19) 

The mode coupling term F2 is given by Equation (45) in 
iBernardeau et all (|2002h . Pre-reconstruction we fix Anl = 
0.19/iMpc"\ which causes damping of the baryonic acous- 
tic feature, and Amc ~ 2.44 which takes into account 
mode coupling. These values are determined by analyzing 
the mean signal of the mocks whilst fixing cz/H/tb and 
Da/vs to the true values ( and n ot using shape parameters). 
See §3 of ISanchez et all (|2013h for a thorough discussion 
of the template and a summary earlier investig ations (e.g, 
ICrocce fc Scoccimarroll2008l : ISanchez et al.ll200y i. 

We compare the results obtained by means of the RPT- 
based model to a popular model denoted as dewiggled, which 
also includes a Gaussian damping of the baryonic acoustic 
feature: 

-Pdcwigglcd(fc, Mfc) = (-PLincar — -PNoWigglc) ' T^{k,fJ.k) 
+ -fNoWiggle; 

(20) 



where the anisotropic damping is defined by: 



V(k, fik 



; exp 



-^e (Misf| + (i-Ml)si) 



(21) 



The PNoWiggie( fc ) is the no-wiggle model given in 
lEisenstein fc Hul ()l998l ) . For a full description of t his model, 



the rea der is referr ed to Eisenstein et al. (|2007bl ) , IXu et al] 



(|2012ah . and §4.3 in lAnderson et al.l(|2013l ). Here we use val- 
ues Ex, Ell = 6, 11 h~^Mpc for the pre-reconstruction case. 



and E 



II 



3h ^Mpc post-reconstruction. 



5.2.2 Post-reconstruction templates 

Equation 1)17^ is used for both the RPT-based and dewiggled 
templates pre-reconstruction. Post-reconstruction templates 
are described in this section. 

Assuming that the reconstruction procedure works 
correctly, one expects, in addition to the sharpening of 
the baryonic acoustic feature, a correction for redshift- 
distortions, yielding an isotropic ^(s). We apply this ap- 
proach in the RPT-based modeling. Due to the sharpening 
of the baryonic acoustic feature, we set /cnl ~ 0.50/iMpc~^, 



which effectively yields the linear ^. The isotropy in the post- 
reconstruction template is introduced by setting cry , /3 = 
and hence ^2 = 0. 

When applying reconstruction we also expect, ideally, 
no need for the coupling term. However, when analyzing the 
mocks, we find that setting ^mc to zero, yields a 0.7% — 1% 
bias in cz/H/r^. For this reason we fix Amc = 2.44 as the 
pre-reconstruction template, which produces lower bias (< 
0.5% see Tables [2l [CTIl . 

For the dewiggled post-reconstruction template we as- 
sume an isotropic Pnl model, but do include ctv = 1 /i~^Mpc 
contributions, which are small at the baryonic acoustic fea- 
ture scale. This i s the same template used in the analysis of 
[Anderson et all (|2013l ). 

As clearly seen in Figure [S] the reconstruction proce- 
dure, as applied on the PTHalos, yields a systematic effect 
in ^2, which reverses sign at scales of the baryonic acoustic 
feature, suggesting there might be an over-compensation of 
the Kaiser effect. We are not concerned by this fact, as we 
are interested in the peak positions to extract cz/H/r^ and 
Da/^s, and not j3. When using each of the templates, linear 
and non linear systematics of the reconstruction procedure 
are corrected by the shape parameters as described in ii5.3.1l 
In i j6.2l we demonstrate that for the RPT-based model the 
post-reconstruction results are essentially unbiassed. 

The RPT-based templates used are plotted in Figures [1] 
(pre-reconstruction) and [2] (post-reconstruction). In the pre- 
reconstruction case we see a clear agreement with the 
and ^0.2. The dewiggle d templates are plotted in Figure 1 of 
[Anderson et~aLl (|2013l) . 



5.3 The model tested 

In this study we focus on the geometric information czjajr^ 
and DA/rs contained in ^ in a model-independent fashion. 
This is done by focusing on the information contained in the 
anisotropic baryonic acoustic feature, and hence marginalize 
over the shape e ffects, i n a similar approach to that used in 
IXu et all (|2012al lbh and [Anderson et al.l (|2012l ). 

For each statistic analyzed we define a model based on 
a template using the following prescription: 



^ mo del / \ 

?stat \Sf) — 0.0 



■^.'^It +Atat(Sf), (22) 



where (^stat = ?||,?±,^o or ^2). The cz/H/rs and Da/ts pa- 
rameters are varied within ^^^^ template application of the 
non-linear AP effect, as described in ij^la'iid Appendix|X] In 
Appendix|B] we also compare the non-linear to the linear AP 
shift and conclude that for DR9-CMASS the linear method 
underestimates constraints by criin<=ar^^non-iinear _ q o 

where is the ID marginalized 68CLr of X=_ff,_DA|!| 



parameters 



5.3.1 The 



As indicated in Equation (|22p . each statistic "stat" is multi- 
plied by its own independent amplitude factor ao stat . These 
factors take into account effective variations of us, galaxy- 
to-matter linear bias, and the effective linear Kaiser boost. 



^ Results from tests on the DR9-CMASS pre-reconstructed ^0,: 
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For each clustering wedge ^n^^ model we add three ad- 
ditional non-linear parameters according to: 



ai 



0.1 



"2 I 

H 1- as II, 

s " 



H ha:; X . 



(23) 



(24) 



When testing for the ^0,2 we apply a similar approach. These 
A{s) terms are applied to the model only after the original 
template is shifted by the AP mapping. Hence, the param- 
eter space used contains ten parameters: 

$10 = [c2/H/r„DA/r„5], (25) 



where 



S — [ao statljfll statl,a2 statljflS statl, 
Ho stat2,ai Stat2,fl2 stat2 , 13 stat2], 



(26) 



where ai statj is the i**^ shape parameter for the 
^—statistic, as described in Equations (|23I24|) . 

In our analysis we find, however, that ao ?2 ii^t well 
constrained both pre- and post-reconstruction (this is not 
the case for the rest of ao stat). We decide to fix this parame- 
ter, and hence are left with a nine parameter space $9, when 
analyzing ^[0,2] • In Appendix[U]we verify that the results ob- 
tained with ^0,2 using <l>g yield similar results (modes and 
uncertainties) to those obtained with ^a/j using $10 both 
pre- and post-reconstruction. In H6.3.3I we describe degen- 
eracies of the shape parameters with cz/H/tb and Da/tb 
constraints. 



5.3.2 Priors 

We limit (DA/rs)/(DA/rs)"'* and ( ffrs)"'^/(Hrs) eac h to 
the region [0.5,1.5]. As suggested bv IXu et al.l (|2012al ). we 
test the effect of applying a Gaussian prior on the warping 
parameter e. We also examine applying a fiat prior on e. 

For most of this analysis, we do not use these priors, but 
we do examine using various e prior values, and report a few 
results with fiat prior |ej< 0.15. These priors are physically 
motivated. First, most reasonable cosmologies would find 
|e|> 0.07 highly improbable. Second, the covariance matrix 
is limited to some extent, and reliability is questionable at 
high dev iations from the fiduc ial cosmology (e.g, see Figures 
16,17 in ISamushia et al.ll2012l '). Third, if the results yield a 
high e, the fiducial cosmology should be revisited. Overall we 
find CMASS results are not sensitive to the choice of prior. 



6 RESULTS 

In this section, we determine the significance with which 
the DR9-CMASS anisotropic baryonic acoustic feature is 
detected, and compare this to simulated realizations. We 
later describe the measurements of cz/H/rB and Dp^/rB. 



6.1 Significance of the detection of the anisotropic 
baryonic acoustic feature 

We generalize the standard technique of determining the sig- 
nificance of the detection of the baryonic acoustic feature to 
the 2D anisotropic case by usage of the clustering wedges. 



and apply this to the DR9-CMASS and the 600 mock real- 
izations. 

The method involves comparing the lowest result of 
a chosen physical mode l to a no-wiggle model . For a no- 
wiggle model we use the lEisenstein fc Hul (|l998h formalism 
(see their §4.2), and derive monopole and quadrupole com- 
ponents using of Equation pop . 

Using this approach as a template, we run the same 
modeling and AP mapping fEauation l22p with the same pa- 
rameter space $10 as the physically motivated templates. In 
the procedure we do not attempt to analyze the clustering 
wedges separately from each other, i.e, we do not attempt 
to quantify significance of detection of the baryonic acoustic 
feature only in or ^x- Instead, we quantify the significance 
of the detection of the anisotropic baryonic acoustic feature 
in the ^(s) by using both ^A/j. This is due to the co- variance 
between the clustering wedges, as well as the strong corre- 
lation between {HrBf'^/{HrB)a.nd (DA/rs)/(L»A/rs)'*'*. AH 
the following results are similar when using the RPT-based 
or the dewiggled template. 

We apply this procedure on both the CMASS and 
the mock catalogs. The results are summarized in Fig- 
ure [S] where the left panels correspond to pre-, and the 
right post-reconstruction. The top two panels correspond 
to the CMASS Ax^ = Xrcf ~ results as a function 
of (HrBf'^/iHrB) and (DA/r,)/(73A/rs)"'^. The thick blue 
lines show the minimum surface of the RPT-based model 
compared to its minimum x?ef • The thin red line corresponds 
to the no-wiggle (no-peak) surface minimum model com- 
pared with Xrcf - The bottom two panels are histogram re- 
sults of the mock realizations, where the CMASS results are 
indicated with the thick vertical line. No priors on e have 
been applied. 

The pre-reconstruction CMASS clustering wedges yield 
a result of (Ax^)min = niin(x?cf ~X^) = 22.2, meaning a 4.7(t 
detection of the anisotropic baryonic acoustic feature, and 
we obtain a similar result after applying reconstruction. This 
result appears to be consistent with the isotropic baryonic 
acoustic feature detection of CMASS-DR9 as reported by 
[Anderson et all (|2012l) . who showed a 5a detection that did 
not improve with reconstruction. 

In the pre-reconstruction case, the CMASS sample ap- 
pears to be on the fortunate side of the mock distribution of 
the detection of the anisotropic baryonic acoustic feature, 
where 68% of the mocks lie between 2.8 — 4.6 a. In the post- 
reconstruction case we see a clear shift of the mocks between 
3.6-5.4 a. 

For later reference we define a subsample of 462 re- 
alisations with a > 3(7 detection as the "> 3a subsam- 
ple", and its complement the "< 3(t" subsample. For a con- 
sistent comparison between the various methods this sub- 
sample is defined when using the pre-reconstruction wedges 
RPT-based method. In the context of the DR9-CMASS vol- 
ume, we find this separation useful for interpretation of the 
(Hr^f'^/iHrB) and {Da/t-b) /{DA/rBf'^ results. For a visual 
of the subsamples in terms of Ax^, please see Figure [S] 

In the following section we analyze how well we expect 
to measure {HrB^'^ /{Hrs) and (DA/rs)/(DA/ra)"'* both 
pre- and post-reconstruction. 
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Figure 5. In the top plots we examine the significance of the detection of the anisotropic baryonic acoustic feature in the CMASS 
clustering wedges by comparing results of two templates: a physical ACDM template (thick blue) and one with no baryonic peak (thin 
red). In each panel in the plots we display the the minimum surface for the marginalized (Hrs)^'^ / {Hrs) (left), and {D\/rs) / [Da/ts)^'^ 
(right). The reference x^ from which each binned result is compared to is that of the best fit of the ACDM model. The left plots correspond 
to the data pre-reconstruction and the right to post-reconstruction. In CMASS we find the significance of the detection of the anisotropic 
baryonic acoustic feature to be ^ Ax^[^ = 4.7 a for both the pre- and post-reconstruction cases. In the bottom plots we run the same 
procedure on 600 mock catalogs and present the histogram of the distribution while indicating the CMASS result. 



6.2 Measuring H, Da- testing methodology on 
mocks 

To test the various assumptions made throughout the 
analysis we first apply the pipeline to our mock catalogs. 
To difTerentiate between systematic effects and peculiarities 
due to mocks with low baryonic acoustic feature signal, in 
Appendix [C] we investigate high S/N mocks to answer the 
following questions (answers based on results in Table [CTjl : 

Does the method outlined in ^5.3\ affect the AP test? 
The RPT-based result entries show that the marginalization 
over the shape information yields small biases (< 0.5%) in 
the geometric information measured. 

Is one ^ template preferred over the other? 
We find that although the RPT-based and dewiggled 
templates yield similar constraints and strong mode corre- 
lations, the dewiggled one yields a ~ 1% bias in measuring 



H^'^/H (Appendix [CT]). The dewiggled template does 
not have a mode coupling term, which might explain 
tendencies to yield more biased mock results than the 
RPT-based template. In ^6.3.3l we report results with varied 
^MC, but defer a more intensi ve investigation of possible 
effects for future studies (e.g. the lTaruva et al]|2010l model). 

Is one ^ combination preferred over the other? 

We find that ^a/j and contain similar constraining power 

(Appendix [C2| . 

Are the resulting distributions of the (Hrs)^'^ /{Hvs) 
and {DA/n)/{DA/rsf'^ Gaussian? 

We find that results of high S/N mocks yield close to Gaus- 
sian (or symmetric) posterior distributions but DR9- volume 
mocks do not. This result is probably due to the fact that 
the DR9 mock volumes contain a large fraction of mocks 
with low S/N anisotropic baryonic acoustic feature. 
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Figure 6. Here we show a quantification of tfie detection of 
the anisotropic baryonic acoustic feature for all 600 mocks, and 
the data both pre- and post-reconstruction. We define the > Scr 
subsample as 462 realizations for which the anisotropic baryonic 
acoustic feature is detected with at least Ax^ = 9 in the pre- 
reconstruction case (grey circles). The complementary are defined 
as the < So- subsample (blue x's). 



Does reconstruction improve/bias the above? 
We find that the reconstructed RPT-based template yields 
a good description of the PTHalo mocks and, on average, 
improves constraints of cz/H/vs and Dj^/rs by ~ 30% 
(Appendix [CH). 

These tests show that the methods applied work well on 
high S/N mocks. Analyzing 600 PTHalo DR9- volumes, we 
find that a non- negligible amount of realizations yield low 
anisotropic baryonic acoustic feature signals. 

Figure [7] shows correlations between H^'^/H and 
Da/D^ modes (top) and their uncertainties (bottom). As 
explained in Appendix [Cl (and apparent in Figure [CTIl . these 
distributions of A{l/H)/{1/ H) and ADa/Da are not Gaus- 
sian. A visual inspection of various individual mocks reveals 
some cases with weak line-of-sight and/or transverse bary- 
onic acoustic features. This is quantified in H6.ll where we 
find that ~ 23% of the realizations have an anisotropic bary- 
onic feature with a significance of less than 3a. For this 
reason, we separate the results to the > 3a subsample (gray 
points) and its complementary < 3a subsample (blue points) 
Note that in both pre- and post-reconstruction the subsam- 
ples are the same as that in pre-reconstruction (for a visual 
see Figure[6|. This separation points to interesting trends in 
H^'^/H and Da/D^ modes and uncertainties. 

Most of the outliers that measure H^'^/H and Da/D^ 
modes at > 10% from the true values tend to be from the 
< 3a subsample in both pre- and post-reconstruction. The 
plot clearly shows that reconstruction substantially improves 
both mode and uncertainty scatters and constraints. 

The uncertainty-uncertainty plots also show that most 
of the extremely large uncertainties are in the < 3a subsam- 
ple. Although post-reconstruction removes the trend differ- 



ences in the uncertainties, we clearly see that the tightest 
constraints are on the > 3a subsample. 

For clarity of the plots and interpretation of results, 
we have applied a |e|< 0.15 prior on the MCMC proposi- 
tions. In the mode-mode plots this limit is shown by the 
dashed lines. The motivation behind this choice is given 
in H5.3.2I Without this prior, we find a systematic "pile- 
up" on the flat prior limit of which is domi- 
nated by the < 3a subsample. We verify that these mocks 
have line-of-sight baryonic acoustic features that are either 
washed out, or contain a with a spurious strong cluster- 
ing measurement at 110 < s < 200 h~^Mpc. For some of the 
"double-mode" realizations (meaning with both at line-of- 
sight baryonic acoustic feature signal and a spurious strong 
feature) the e prior strengthens the true mode. For realiza- 
tions with strong spurious features the e prior causes them to 
move from H*^'^/H = 0.5 closer to the e — —0.15 boundary. 

All the above trends appear in both templates examined 
(RPT-based, dewiggled), and in both clustering wedges and 
multipoles. 

In Table [2] we summarize the mock results of H^'^/H 
and Da/Da^ modes and uncertainties and their scatter. 
Most entries are for the RPT-based clustering wedges pre- 
and post-reconstruction. For completeness, the first and last 
entries include the dewiggled templates as well as including 
results of multipoles in all templates. The sample examined 
is indicated (e.g., full sample or the > 3a subsample) as well 
as if a prior on e is used. For example, we investigate vari- 
ous jej priors, or restricting to realizations with H^'^/H and 
Da/Da^ modes within 14% from the true values, or both. 

Regarding the post- reconstruction RPT-based ^||,x we 
notice that H^'^/H has in all cases a median mode bias of 

< 0.3%, and Da/D%'^ < 0.1%. Pre-reconstruction mode re- 
sults, on the other hand, improve substantially when apply- 
ing the various priors and cuts {> 3a subsample, |e|< 0.15, 
mode limitation). These results show the effects of mocks 
with low anisotropic baryonic acoustic feature signal. For ex- 
ample when limiting the sample to the most constrained 2/3 
of the realizations (meaning 394/600), the bias on H^'^/H 
improves from 3% to 0.4%, and of Da/D%^ from 0.6% to 

< 0.1%. 

The H^'^/H and Da/ uncertainties improve in dif- 
ferent manners when applying the various priors and cuts. 
The most noticeable trend, which is common for both pa- 
rameter results, is the reduction of the scatter on the uncer- 
tainty when applying the |e|< 0.15 prior. 

For ill-constrained DR9 volumes the median uncertain- 
ties vary with choice of e. On the other hand, for well- 
constrained realizations, such as CMASS-DR9, results do 
not depend on the e prior (see i)6.3p . 

We also find that the dewiggled pre-reconstruction tem- 
plate yields similar H^'^ / H and Da/D^ constraints as the 
RPT-based ones, although the dewiggled pre-reconstruction 
template shows a systematic bias of ~ 1% on Da/D^ . This 
effect is not apparent in the high S/N mocks (Appendix lC2|) . 
which yield a median 1.4% bias on H^'^/H, which is not ap- 
parent here. The post-reconstruction dewiggled wedges re- 
sults are in line with the RPT-based. 

Perhaps the most notable feature in Table [2] is that the 
scatter in the H^'^/H modes is different from the median 
of the uncertainties. Focusing on the most constrained sub- 
sample (the bottom entry), we see that the scatter in the 
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Figure 7. Pre-(Left) and post-reconstruction (Right) distributions of Q|| = (Hrs)^'^ / (Hrs) and ctx = {Da/ts) / (Dj^/ts)^'^ modes and 
their uncertainties of the mock PTHalos using the RPT-based The top panels show the scatter of mode measurements; the 

bottom presents the scatter of uncertainties. Each panel presents the results of all 600 mock realizations, where the grey dots are the 
> 3(T subsample (462 realizations), and blue for the complementary < So" subsample. The solid contours in each panel are the 68, 95% CL 
regions for the > Scr subsample (gray) and the full sample (black). The cross-correlation coefficient for the > 3cr subsample in each panel 
is indicated by r. Numerical results are summarized in Table |2] In the top panels we emphasize the constant a and e lines, as indicated 
(where the thicker line of each indicates the larger value). In the bottom panels we mark the DR9-CMASS uncertainty measurement 
(red filled squares). For plotting purposes we apply a prior of |e|< 0.15. 



H^"^ /H modes is smaller than the median of the uncertain- 
ties in all cases. In ii6.4l and Appendix [Cl we show that this 
should improve with higher S/N samples. For Da/ we 
see that the scatter of the modes and median of the uncer- 
tainties are fairly similar. 

The fiducial cosmology of the analyses is the true cos- 



mology of the mocks. We defer testing possible effects of 
using an inco rrect fiducial cosm ology (for preliminary tests 
on mocks see iKazin et al.]|2012h . 

To summarize, we find that a significant minority of 
DR9-CMASS pre-reconstruction realizations yield unreli- 
able results. However, the majority > 3a subsample yields a 
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Table 2. Mock DR9 PTHalo results 



5 (# of realizations) 



Aa||/a| 



Aax/c«x 



Full sample, no priors: 

RPT-based wedges pre-Rec (600) 0.970 ±0.188 0.132 ±0.075 1.006 ±0.087 0.048 ± 0.061 

RPT-based wedges post-Rec (600) 0.997 ± 0.101 0.068 ± 0.065 1.000 ± 0.042 0.034 ± 0.033 

RPT-based multipoles pre-Rec (600) 0.986 ±0.194 0.102 ±0.076 1.001 ± 0.082 0.050 ± 0.053 

RPT-based multipoles post-Rec (600) 0.992 ±0.176 0.077 ±0.083 1.002 ±0.052 0.037 ± 0.024 

dewiggled wedges pre-Rec (600) 0.983 ± 0.190 0.129 ± 0.075 1.014 ± 0.086 0.047 ± 0.060 

dewiggled wedges post-Rec (600) 0.999 ±0.106 0.065 ± 0.067 1.002 ±0.047 0.033 ± 0.036 

dewiggled multipoles pre-Rec (600) 0.990 ± 0.183 0.100 ± 0.072 1.008 ± 0.086 0.049 ± 0.047 

dewiggled multipoles post-Rec (600) 1.002 ± 0.134 0.056 ± 0.077 1.000 ± 0.045 0.030 ± 0.025 

Full sample, |e < 0.15|: 

RPT-based wedges pre-Rec (600) 0.982 ± 0.114 0.098 ± 0.049 1.002 ±0.050 0.044 ± 0.034 

RPT-based wedges post-Rec (600) 0.998 ± 0.067 0.064 ± 0.038 0.999 ± 0.038 0.033 ± 0.020 

> 3fT subsample, no priors: 

RPT-based wedges pre-Rec (462) 0.983 ±0.146 0.103 ±0.065 1.003 ±0.064 0.042 ± 0.044 

RPT-based wedges post-Rec (462) 0.997 ± 0.086 0.061 ± 0.061 1.000 ± 0.034 0.032 ± 0.026 

> 4.0(7 subsample, no priors: 

RPT-based wedges pre-Rec (208) 0.990 ±0.111 0.074 ± 0.054 1.001 ± 0.043 0.034 ±0.027 

RPT-based wedges post-Rec (208) 0.996 ± 0.079 0.053 ± 0.051 0.999 ± 0.030 0.028 ± 0.020 

> 4.5cr subsample, no priors: 

RPT-based wedges pre-Rec (104) 0.992 ± 0.099 0.065 ± 0.050 1.000 ±0.035 0.031 ± 0.021 

RPT-based wedges post-Rec (104) 0.999 ± 0.045 0.052 ± 0.040 0.997 ± 0.024 0.028 ± 0.009 

> 3(T subsample, |e < 0.15|: 

RPT-based wedges pre-Rec (462) 0.988 ±0.102 0.089 ± 0.042 1.001 ± 0.048 0.040 ± 0.023 

RPT-based wedges post-Rec (462) 0.997 ± 0.061 0.059 ± 0.035 0.999 ± 0.031 0.031 ± 0.014 

> 3o- subsample, |e < 0.15|, |1 - mode|< 0.14: 

RPT-based wedges pre-Rec (394) 0.996 ± 0.060 0.087 ±0.035 1.000 ±0.037 0.040 ± 0.021 

RPT-based wedges post-Rec (450) 0.998 ± 0.046 0.059 ± 0.032 0.999 ± 0.029 0.031 ± 0.014 

RPT-based multipoles pre-Rec (374) 0.999 ± 0.061 0.079 ± 0.030 0.997 ±0.038 0.044 ±0.011 

RPT-based multipoles post-Rec (434) 1.001 ± 0.047 0.062 ± 0.030 0.998 ± 0.031 0.033 ± 0.011 

dewiggled wedges pre-Rec (392) 1.003 ± 0.063 0.087 ± 0.036 1.009 ± 0.038 0.039 ± 0.019 

dewiggled wedges post-Rec (445) 1.003 ± 0.047 0.056 ± 0.034 1.000 ±0.031 0.030 ± 0.016 

dewiggled multipoles pre-Rec (371) 1.003 ± 0.061 0.077 ± 0.030 1.006 ± 0.037 0.042 ± 0.019 

dewiggled multipoles post-Rec (434) 1.006 ± 0.043 0.051 ± 0.033 0.999 ± 0.028 0.028 ± 0.007 

* The Oji and ox columns show the median and rms of the modes. 

* The A«|[/a|| and Aax/^x columns show the median and rms of the fractional uncertainties. 



low bias result (< 0.5%). Moreover, the results show that the 
post-reconstruction wedges results yield low bias (< 0.3%) 
with both RPT-based and dewiggled. We also find that for 
a DR9 volume we expect non-Gaussian likelihood profiles of 
cz/H/ts and Da/j's in both pre- and post-reconstruction. 
We next turn to apply the same method used here on the 
data both pre- and post-reconstruction. 



accuracy). The correlation coefficient betwen cz/H/rs and 
Da /''a is measur e d at —0.5, similar to that predicted by 
ISeo &i Eisensteiiil (|2007h . The best fit model shows an ex- 
cellent fit at x^/dof=0.82 with dof=66 degrees of freedom. 
Compared to the mocks realizations, this result is better 
than 398/600 mocks. With the pre-reconstruction we 
obtain /'^oi—Q.GA (better than 578/600 realizations). 



6.3 DR9-CMASS H, Da results 

In this section we present our measurements of cz/H/rs and 
Da/ts in the DR9 CMASS dataset. 

Our main pre- and post-reconstruction results are sum- 
marized in Figures [2] and |8l Post-reconstruction we mea- 
sure cz/H/ra = 12.28 ± 0.82 (6.7% accuracy; uncertainties 
are quoted at 68% CL) and Da/^s = 9.05 ± 0.27 (3.0% 



Figure |8] compares the posterior results (solid red lines) 
with a Gaussian approximation (dashed blue lines) , based on 
the same quoted modes, uncertainties and cross-correlation 
coefficients. In both the pre- and post-reconstruction cases, 
we see that the Gaussian approximation describes the 
68.27% CL region fairly well, but clearly underestimates 
the 95.45% CL reigon. We also note that the full pos- 
terior 99.73% CL regions obtained both pre- and post- 
reconstruction are not well defined. These indicate the lim- 
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Figure 8. CMASS results prc-roconstruction (left) and post (right). The marginalized results of cz/H/rg (right panels) and Dj^/r^ 
(top panels), and the joint constraints (bottom panels). The solid red lines are the posterior, and the dashed blue lines are a Gaussian 
approximation, as described in the text. The panels indicate the modes, 68% CL region boundaries, proposition-mean, proposition 
standard deviation, skewness and cross-correlation coefficient (r). The contours indicate the 68.27,95.45% CL regions. For plotting 
purposes the post-reconstruction likelihoods assume a prior |e|< 0.15. The gray dashed lines indicate the fiducial cosmology. 



ited S/N in these measurement. We expect the agreement 
to improve with larger samples. 

In the top plot of Figure [9] we show a direct compar- 
ison of the likelihood profiles pre- and post-reconstruction 
of the RPT-based clustering wedges. Both results appear 
to be similar, well within the 68% CL region, although in 
the the post-reconstruction case cz/H/rs is not as tightly 
constrained. 

For an average mock DR9- volume realization we find a 
mode cross-correlation of ri/[j,^Dj^ ~ 0.35 — 0.4 (or ~ 0.6 
when examining the high S/N mocks; Appendix ICSP should 
be expected, where ri/u is the cross-correlation between 
the cz/H/rs modes obtained when using one method (here 
pre-reconstruction) and when using a second (here post- 
reconstruction), and similar for r_D^, when discussing Da/ts 
results. Also, although one does expect tighter constraints 
when applying reconstruction, the DR9 mocks indicate a 
19% (116/600) possibility of not improving cz/H/vs- Using 
mocks with expected S/N of the final BOSS footprint (de- 
scribed in H6.4p . this probability is reduced to ~ 1.5%. 

The CMASS cz/H/vs, Da/ts results are summarized in 
Table |3] along with various related parameters. 

6.3.1 Comparing results of various ^ methods 

The results quoted in the previous section are obtained when 
using the A/i = 0.5 clustering wedges with the RPT-based 
template. Table [3] contains the results obtained for eight dif- 
ferent combinations of statistics. 

When applying the dewiggled template we obtain simi- 
lar results to those obtained with RPT-based one. According 
to our mocks we expect r^iH^TDA ~ 0.5 — 0.65 amongst the 
templates both pre- and post-reconstruction. 

We apply the same test on the [^0,^2] multipoles and 
obtain slightly different results, but consistent within the 
68% CL regions, as seen in the bottom plot of Figure [§1 Ac- 
cording to the DR9 mock realizations we expect cross corre- 



lations between wedges results to multipoles by ri/£f,r_D^ ~ 
0.4-0.45. 

Figure [TD] displays cz/H/rs, DA/rs likelihood profiles 
of all eight different methods analyzed here. The plot shows 
that all methods yield consistent results. The ^0.2 pre-rec 
(both RPT-based and dewiggled) cz/H/rs profiles appear to 
be wider than the rest, where the ^0,2 post-rec (both RPT- 
based and dewiggled) appear to be the furthest from the rest, 
although clearly consistent within the 68 — 95% CL regions. 
These differences are as expected based on the results from 
the mocks (for a visual of higher S/N mock results see top 
plot in Figure [12)) . We investigate various methods of shape 
parameters, and find similar results. 



6.3.2 Robustness of results to the range of fitted scales 

As discussed in tj5.3l these measurements focus on the infor- 
mation of the anisotropic baryonic acoustic feature and not 
from the full shape. As such, we do not expect dependency 
of our results on the range of scales used in the analysis. 

The results quoted in the previous sections are obtained 
when analyzing data in the region of separations between 
[smin, Smax] ~ [50, 200]. We Compare the results obtained for 
various choices of Smin, Smax. Figure[TT]shows the comparison 
of the results. 

We find that, for the most part, the range of analy- 
sis does not affect our main results: mode values, uncer- 
tainties, cross-correlation coefficient or skewness. Regions of 
exception involve those with Smin > 65/i~^Mpc, in which 
the cz/H/ra uncertainties increase from ~ 6% to 7% and 
even higher, when limiting to Sniax=160 fe~^Mpc. This re- 
sult could be explained by the fact that in this latter test 
the full dip of the baryonic acoustic feature is not used, 
and shape parameter values that cause spurious dips are 
accepted, whereas for lower values of Smin they are not. 
We conclude that a more reliable result would include data 
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Table 3. CMASS DR9 (z) = 0.57 results 



t, 


a. 


CZ f Tl f Vs 








TI 

£1 




No prior on e: 














km-s^^Mpc"! 


Mpc 


RPT- based 5|| ^ pre-Rec 


1.042 


12.41±0.75 (6.1%) 


1.006 


8.92±0.27 


(3.0%) 


-0.50 


89.9±5.6 


1367±45 


RPT-based ^0.2 pre-Rec 


1.072 


12.77±1.15 (9.0%) 


0.989 


8.77±0.36 


(4.1%) 


-0.72 


87.4±7.9 


1344±57 


dewig C||,x pre-Rec 


1.055 


12.57±0.73 (5.8%) 


1.014 


8.99±0.28 


(3.1%) 


-0.57 


88.8±5.3 


1378±46 


dewig ^0,2 pre-Rec 


1.070 


12.74±1.06 (8.3%) 


1.008 


8.94±0.33 


(3.7%) 


-0.72 


87.5±7.4 


1370±53 


RPT-based 5|(,x post-Rec 


1.031 


12.28±0.83 (6.8%) 


1.020 


9.05±0.25 


(2.8%) 


-0.51 


90.8±6.2 


1386±42 


RPT-based ^o,2 post-Rec 


0.974 


11.60±1.44 (12.4%) 


1.055 


9.36±0.34 


(3.6%) 


-0.67 


96.2±12.0 


1434±55 


dewig 5||_x post-Rec 


1.026 


12.22±0.96 (7.8%) 


1.020 


9.05±0.24 


(2.7%) 


-0.54 


91.3±7.3 


1386±41 


dewig ^0,2 post-Rec 


0.974 


11.60±0.79 (6.8%) 


1.046 


9.28±0.27 


(3.0%) 


-0.62 


96.2±6.7 


1422±45 


prior |e|< 15%: 


RPT-based ^\\^± pre-Rec 


1.042 


12.41±0.75 (6.1%) 


1.006 


8.92±0.27 


(3.0%) 


-0.50 


89.9±5.6 


1367±45 


RPT-based ^0,2 pre-Rec 


1.072 


12.77±1.15 (9.0%) 


0.989 


8.77±0.36 


(4.1%) 


-0.75 


87.4±7.9 


1344±57 


dewig ?||^x pre-Rec 


1.055 


12.57±0.72 (5.8%) 


1.014 


8.99±0.27 


(3.0%) 


-0.53 


88.8±5.2 


1378±45 


dewig ^0,2 pre-Rec 


1.070 


12.74±1.06 (8.3%) 


1.008 


8.94±0.33 


(3.7%) 


-0.72 


87.5±7.4 


1370±53 


RPT-based 5||,x post-Rec 


1.031 


12.28±0.82 (6.7%) 


1.020 


9.05±0.27 


(3.0%) 


-0.50 


90.8±6.2 


1386±45 


RPT-based ^0,2 post-Rec 


0.974 


11.60±0.87 (7.5%) 


1.052 


9.33±0.36 


(3.8%) 


-0.78 


96.2±7.3 


1430±57 


dewig 5||jX post-Rec 


1.026 


12.22±0.91 (7.4%) 


1.020 


9.05±0.27 


(3.0%) 


-0.53 


91.3±6.9 


1386±45 


dewig ^0.2 post-Rec 


0.974 


11.60±0.73 (6.3%) 


1.046 


9.28±0.33 


(3.6%) 


-0.63 


96.2±6.2 


1422±54 



* We define a|| = (Hr^)^"^ / {Hr^) and ax = (Da As)/(DaAs)*^'*. 

* Uncertainties A quoted correspond to half of the 68% marginalized CL region (68CLr). In parentheses is the mean percentage of 
68CLr. 

* All values are unitless, unless otherwise indicated. 

* Fiducial values used at (z) = 0.57: {cz/H/r^Y = 11.93, (Da/t-s)' = 8.88, based on WMAP5 co smology jKomatsu et al.ll2009l ). 

* The iT(0.57), Da (0.57) columns assume WMAP5 result: rs(zd) = 153.3 ± 2.0 Mpc (Table 3 in lKomatsu et al.ll2009i v 

* ^a||,a^ is the cross-correlation coefficient for cz/H/ts and Da/tb- 



points along the full shape, even though that information is 
marginalized over through the linear bias and A{s) terms. 

We do not consider analyses with Smin < 50/i~^Mpc, 
because the templates used do not describe well the velocity- 
dispersion damping in the PTHalo mock-mean signal, and 
hence models would too heavily depend on the A{s) terms. 

In all ranges investigated the x^/dof is between 0.6—0.8, 
with the Smax = 180 h~^Mpc yielding the best fits, although 
not significantly better ones. 

6.3.3 Regarding the nuisance and fixed parameters 

As described in i\5.3l we use a set of ten parameters $10. To 
best understand the effects and correlations of these parame- 
ters amongst themselves and with cz/H/rs, Da/j's we exam- 
ine the results of both the data and the mock-mean signal. 
We perform these tests both pre- and post-reconstruction in 
both templates for and ^0,2. 

Overall, we do not see particular strong correlations 
between the A{s) shape parameters with cz/H/ra, Da/ts, 
where most cross-correlations are r < 0.2, but do illuminate 
a few findings of interest. 

Most of the shape parameters have marginalized likeli- 
hood profiles that are fairly symmetric (low skewness). We 
find that amplitude parameters Oq || and ± are uncorre- 
lated with each other. All correlations of these parameters 
with cz/H/va and DA/va are r < 10%. The constant param- 
eters (03) are uncorrelated to cz/H/ra and Da/j's, as ex- 



pected. The other shape terms have weak correlations with 
cz/H/ra and Da/»-s, (at r < 0.2). 

The most important finding of the shape parame- 
ters, however, regards the ao (2 (the amplitude of the 
quadrupole). In both pre- and post-reconstruction its 
marginalized likelihood profile is not well constrained, caus- 
ing strong skewness in the joint likelihoods with other pa- 
rameters. We decide to fix its value, which yields results 
similar to ^A/j, where this behaviour is not present. 

We find all the above similar for the data and 
mock-mean in the pre-reconstruction case. In the post- 
reconstruction case this is true as well, after we apply a prior 
|e|< 0.15. Before applying the prior, the 99.7% CL region is 
not well defined as the MCMC chains tend to accept values 
at the low limit set (Hra)^"^ /{Hra) = 0.5. 

Finally we address the question of the Amc pa- 
rameter in the RP T -base d template (Equation llSp . 
ICrocce fc Scoccimarrd (|2008l ) introduced this parameteri- 
zation to effectively take into account the coupling be- 
tween the fc— modes, which results in a 0.5% shift in the 
peak position in ^o- To obtain reliable templates of the 
post-reconstruction ^a/j and we find that a model with- 
out an Amc term yields biassed results in the mocks, by 
about ~ 1% in {Hra)^'^/{HrB). When analyzing the post- 
reconstruction CMASS ^[|_x results, we see a shift in a from 
1.026 {Amc = 2.44) to 1.030 (Amc = 0), a 0.4% increase. 
The 1 + e value is similar at 1.003. This results in a 0.3% 
shift in cz/H/ra and 0.2% shift in Da/ts, well below the un- 
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Figure 9. Comparison of tiie CMASS cz/H/ts and Dp^/r^ re- 
sults obtained witii the pre-reconstruction wedges witli alterna- 
tive methods. The top plot shows a comparison with the post- 
reconstruction result. The Bottom plot shows a comparison 
with the pre-reconstruction clustering multipoles ^0i€2- AH meth- 
ods use the RPT-based template. The contour plots show the 
68, 95% CL regions. The solid lines are the fiducial cosmology. 



certainties. For the post-reconstruction ^0.2 we find similar 
results. 



6.4 Final CMASS Forecasts 

By the conclusion of BOSS (2014), the survey will cover 
three times the area of the data set analyzed here, meaning 
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Figure 10. Comparison of the CMASS-DR9 cz/H/rs and Da/t-s 
marginalized profiles obtained with all the methods tested here. 
The top panel shows results when using the and the bottom 
panel when using l;o,2- The contour plots show the 68% CL re- 
gions. The solid lines are the fiducial cosmology. To guide the eye 
wc plot the regions of constant a and e, as indicated in the legend 
(where the thicker line of each indicates the larger value). 



The full CMASS sample will have a volume three times as 
large. By stacking the PTHalo mocks by groups of three, we 
can effectively, to first order, forecast the cz/H/rs, Da/ts 
results of the full CMASS galaxy sample. Using the 600 real- 
izations, we analyze here results of 200 ^||,x stacked mocks. 

It is important to emphasize that the estimates yielded 
here should be considered maximum bounds. We argue this 
due to the fact that the dj used is the same DR9 vol- 
ume covariance matrix as in Equation (|15|l but divided by 
three. This means that we do not account for noisy cross- 
correlations which should be reduced with the actual full 
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Figure 11. This plot shows that the CMASS RPT-based g|| _l 
pre-reconstruction results are insensitive to the range of analysis 
used when Smin < 65. The x axes are the minimum separation 
used Smini where we compare results of Sj^jn =50,55,60,65,70,75 
h~^Mpc with maximum separations of Sniax =160 (black circles), 
180 (red crosses), and 200 /i~^Mpc (blue squares). The number of 
degrees of freedom (dof) and x^/dof are quoted. All uncertainties 
indicate the 68% CL regions. The blue dot-dashed lines are the 
fiducial input values used to convert z into comoving distances, 
and the black dashed lines are the chosen result quoted in Table 
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CMASS geometry, thus we expect the constraining power 
to be tighter when using a more reliable dj. Furthermore, 
we note that replicating the DR9 geometry does not improve 
the reconstruction boundary effects. 

Figure [121 displays the cz/H/rs and Da/vs results ob- 
tained by means of the expected modes and uncertainties, 
comparing between post- and pre-reconstruction ^||_x (top), 
and post-reconstruction ^||_^ to ^0,2 (bottom). 

When comparing cz/H/vs, Da/tb results of the ^o_2 
to the ^||,x we find strong correlations where mode bi- 
ases are sub 0.3%. Uncertainties show that no method is 
preferred over the other. When comparing pre- and post- 
reconstruction wedges, we find a r ~ 0.52 between the 
modes. 

When applying reconstruction, the cz/H/rs uncertain- 
ties are predicted to improve from 0.045 ± 0.017 to 0.030 ± 
0.006, a 33% improvement. For Da/tb the improvement is 
forecast to be from 0.024 ± 0.007 to 0.017 ± 0.003, a ~ 30% 
improvement. The mock result distributions yield Gaussian- 
like features, although application the K-S tests indicates 
they are not Gaussian. These trends are similar to those 
seen with 100 six-stacked mocks (see Appendix [H)) . 



7 DISCUSSION 

The cz/H /tb and Da/ts results obtained here are consistent 
across the various techniques investigated: 

(i) Cll.x, Co,2 

(ii) ^ template: RPT-based, dewiggled 

(iii) pre- and post-reconstruction 
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Figure 12. a,, = (HrBf-^ /{Htb) and aj_ = {DA/rB)/{DA/rBf'^ 
mode and fractional uncertainty forecasts of 200 pscudo final 
BOSS CMASS volumes. In all plots the y-axis results are for 
post-reconstruction wedges. In the top plot the x-axis results are 
for pre-reconstruction wedges, on the bottom post-reconstruction 
multipoles. In each plot the comparisons are between q:|| modes 
(top left panels), a± modes (top right), Aa||/o|| uncertainties 
(bottom left), Aa±/aj_ uncertainties (bottom right). The cross- 
correlation in each is r. The dashed red lines are the 68% CL 
regions. For the comparison, the red boxes are the DR9-CMASS 
results. 
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The likelihood profiles obtained with these eight combina- 
tions investigated, are shown in Figure[lO] (as well as Figure 
[9] and Table [3]). Differences between the results are as ex- 
pected from mock simulations. 

As these posteriors are not Gaussian, we provide joint 
2D marginalized likelihood profiles of cz/H/rs and Da/ts, 
as well as provide the CMASS-DR9 ^||_x and dj, C,^^ on 
the World Wide Web0 We conclude this study by using 
results obtained post-reconstruction over those yielded pre- 
reconstruction, because we show that mock results expect 
an improvement of 30% in the marginalized constraints of 
cz/H/va and Da/vb, even though this is not the case in 
the data. We also prefer the RPT-based template over the 
dewiggled due to the larger bias in the mock results when 
using the latter. In the data, we find the posteriors to be 
similar regardless of choice of template (see Figure [TO)) . 

Comparison of our results to other analyses of the 
same data set can be found in the following studies. 
[Anderson et al.l (|2013l ) measures cz/H /r^ and Da/ts by ap- 
plying a similar model-independent method on the ^0,2, us- 
ing the same dewiggled templates. The main differences in 
analysis involve their use of a grid of a and e, where the 
rest of the nuisance parameters are determined by the least- 
squares method. We perform extensive comparisons between 
the methods, and find the cz/H/vs a nd Da/tb results to b e 
fairly similar (se e Figu res 13 and 14 in lAnderson et al .l2013h . 
[Anderson et all (|2013l ) continue to use results obtained in 
both studies to produce a "consensus result" and calculate 
cosmological implications. 

Model-depe ndent analyses are performed on the full 
shap e of f || 1 (jSanchez et al.l 20131') and ^0,2 (|Reid et all 
l2012l : IChuang et al.ll2013l '). [Sanchez et all (|2013l ') shows that 
results among s t thes e studies are compatible. Figure 15 in 
ISanchez et aD (HoTl) shows a comparison between our pre- 
reconstruction model independent result and their results 
from the full shape which are independent of parameter 
space, but assume / follows GR predictions. They find an 
excellent agreement with our results, although tighter con- 
straints as the baryonic acoustic feature only method effec- 
tively accepts parameter values (e.g, Q,m) that the full shape 
does not. 



The application of reconstruction leads to a significant im- 
provement of detection of the peak in mock catalogues from 
3.7(T ± 0.9(7 to 4.5(T ± 0.9(7 (median ± standard deviations; 
see Figure [SJ. Pre-reconstruction mocks also show that 23% 
(138/600) yield a detection lower than 3(7, whereas post- 
reconstruction 4.6% do (28/600; FigurefSjl. Although we see 
clear improvement in the average mock realization, the sig- 
nificance of the detection of the anisotropic baryonic acous- 
tic feature in the data does not improve after applying re- 
construction. We find this, however, consistent with 89/600 
(15%) of the mock realizations (Figure [6l). 

To obtain geometrical constraints that are model 
independent, we use information from the post- 
reconstruction anisotropic baryonic acoustic feature 
and measure cz/H/rs = 12.28 ± 0.82 (6.7% accuracy) and 
Da/ts = 9.05 ± 0.27 (3.0%) with a correlation coefficient 
of —0.5 (uncertainties are quoted at 68% CL). In terms of 
constraining cz/H/tb and Da/tb, the pre-reconstruction 
DR9-CMASS yields mutual constraints tighter than 584/600 
of the mocks, putting it in the fortunate top 2.5%. In the 
post-reconstruction case this is reduced to the top 444/600, 
meaning the top 26%. Ahhough CMASS-DR9 results do 
not improve with reconstruction, mock catalogs indicate 
that, on average, one should expect an improvement of con- 
straining power of ~ 30%. Throughout this study we show 
that the posteriors of cz/H/tb and Da/tb from the DR9 
volume are not expected to be Gaussian. In i|7] we explain 
how to use the results presented here, pointing out that the 
provided full likelihood fun ction should be used i nstead of 
a Gaussian approximation. [Anderson et al.l (|2013l ') analyze 
cosmological consequences of this measurement. 

In our analysis of mock catalogues we also demonstrate 
that the constraining power of ^0,2 and i\\,_L are expected 
to be similar. With this information we conclude that the 
analysis of the clustering wedges and comparison to the mul- 
tipoles technique, as performed here, is vital for testing sys- 
tematics when measuring cz/H/tb and Da/tb- Here we use 
wide clustering wedges of A/i = 0.5, which are fairly cor- 
related (see Figure [3} . As long as covariances can be ade- 
quately taken into account, this method could be general- 
ized to narrower clustering wedges, as future surveys will 
yield better signal-to-noise ratio. 



8 SUMMARY 

In this study we investigate the ability of the BOSS DR9- 
CMASS volume to constrain cosmic geometry at 2 = 
0.57, through the use of the AP technique applied on the 
anisotropic baryonic acoustic feature. We analyze the infor- 
mation contained in the anisotropic baryonic acoustic fea- 
ture, for the first time, using a new technique called clus- 
tering wedges ^A/j, and compare results to the multipoles 
Co, 2. 

We find the anisotropic baryonic acoustic feature to 
be detected in DR9-CMASS at a significance of 4.7(7 com- 
pared to a featureless model fi ]6.1|) . We find this level to be 
fairly fortunate (from a cosmological variance perspective), 
but consistent with that expected from mock realizations. 



|http: //www. sdss3 ■ org/science/boss_publlcations .php 
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APPENDIX A: THE AP MAPPING IN 
PRACTICE 

Here we describe the geometrical correction mapping (or AP 
shifting) of ID statistics as ^||,x,Co,2- 

As we compare a ^ template to data which are affected 
by geometrical distortions we must distinguish between two 
sets of coordinate systems, which are, ultimately, related 
through H and Da- 

In 31] we define the geometric distortions of the compo- 
nents of s. In the final product, though, we use its absolute 
value and fi related by: 



(Al) 



where S|| is the line-of-sight separation component. 

The template, from which the model is constructed, 
is calculated in a "true" or "test" coordinate system St , /xt , 
where the data are in shifted axes based on the fiducial cos- 
mology, hence we define its separations and angles Sf,^f. 
Because we apply the model to the data, the model, which 
is based on the template, should be in the fiducial coordi- 
nate system, as well, hence the AP shifting of the template 

to C*"""'''"*'=(sf,Mf)- 

Using Equations (IS])-® along with Equations (|A1|I we 
obtain 



St = Sf ■ y^Q||M? + ai(l-Mf), 



and 



^J.t = /if 



y'^iiMf + aili-M?) 



(A2) 



(A3) 



After ^'°™P'''"'(sf,/if) is produced (see below for details of 
its construction), we calculate: 

^template / \ J 

^ (si,iii)dfii 



template/ \ 

Cam = 



(A4) 



for the clustering wedges. For the niultipoles we calculate: 



template 



(sf) = (2£+l 



^0 



emplate / \ /-• / \ t 

(A5) 

To calculate ^"'"'^''''''(sf > Mf) in practice we apply the 
following steps: 

(i) At every point of the integration we use Equations 
(|A2[) - (|A3|) to convert the fiducial Sf,/if into the template 
true coordinates St,Mt- 

(ii) We interpolate stored arrays of a pre-calculated Co, C2 
templates to the resulting St value. For details regarding the 
templates used see tj5.2l 

(iii) We calculate C(sf,/if) by interpolation of 
C (st (sf,Q||,ax,^tf) ,/it (a||,Qx,/if))=Co(st) + /l2(Mt)C2(st) 

Note that to calculate template ^gqyg^^jQjj |A5p we 
need to calculate £.2{fJ.i), where for the f Eg nation IA4p 
this is not needed. We test our algorithm by applying it on 
mock catalogues in which we assume an incorrect fiducial 



cosmology, and apply the above algorithm and obtain the 
true 1/H and Da values. 

In this method we make two main assumptions. First, 
the AP shifting is based on a template that consists of mul- 
tipoles £ — 0,2. This template can be easily expanded to 
higher orders of £, although at scales of interest ^ > 4 com- 
ponents should be fairly weak. Second , we assume th e plane 
parallel approximation for each pair. IWagner et all (|2008l ) 
sh ow that light - conin g yields minimal effects at « = 1, 3, as 
do lKazin etal] (|2012l ') at z = 0.35. 



APPENDIX B: LINEAR VS. NON-LINEAR AP 
EFFECT 

Throughout this analysis we apply the non-linear AP cor- 
rection as described in Appendix [X] In this section we in- 
vestigate differen ces with the linear AP effect as used in 
Xu et all (l2012b|). This linear approach was introduced in 
Padmanabhan fc White! (12008'') in the P{k) formulation, and 
analyz ed in ^ in iKazin c t al. ( 2 01^. However, as pointed 
out bv lPadmanabhan fc White! (|2008l ). this linear approach 
breaks down when |e|> 2%, which is clearly the case in the 
DR9-CMASS for a large part of the 95% CL region. 

The linear AP correction, when applied on the cluster- 
ing multipoles, is as follows: 



?o(st) 



Co(«Sf) + 

'2^^2(2; 



^ / x — asf 



7 J 7 dln{x) 



+ 



2e 



d^o{x) 



dln(x) 



(Bl) 



(B2) 



Here we neglect terms of order 0{e^), as well as ^4 
terms. (For a discu ssion of higher order terms see §2.2.4 in 
iKirkbv et al]|2013l .') 

The left plot of Figure IBll shows the results obtained 
when applying the non-linear AP (thick blue) and the linear 
correction (thin red) as to the CMASS-DR9 The dotted 
and dashed lines convey constant values of a and e, respec- 
tively. 

The results clearly show that the linear correction 
under-estimates the uncertainties of cz/H/va and Da/tb 

by crl^"<=='7a^°"-""'==''- = 7.2/9.6 and f^l-ear/^non-linear _ 



3.2/3.9, where al" 



'-Da - 
is the 68CLr of X=H,Da. The 
method results agree fairly well where e is small (and regard- 
less of a) , but differ as e grows. These differences should vary 
with the choice of the fiducial model, as well as the volume 
investigated. 

We apply a similar comparison for a mock-mean signal 
(of 600 mocks) with the dj divided by three (as in H6.4p and 
plot the results in the right of Figure [Bl] In this higher S/N 
test we clearly see that the two methods agree with each 
other extremely well, due to the fact that e is low. There 
is a slight under-estimation of the linear approximation at 
the 95% CL region. Note that here we test the case where 
the fiducial H and Da correspond to the mock true values 
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Figure Bl. We use the RPT-based multipoles pre-reconstruction to test the hnear (thin red) AP correction against the non-linear 
(thick blue) in constraining cz/H/tb and Dx/ts- The left plot shows results for CMASS DR9 investigated here (contours are 68, 95% CL 
regions), and the right plot for projections of the final BOSS CMASS footprint (contours are 68,95,99.7% CL regions). (As mentiond 
in i]6.4l this forecast should be considered overestimated constraints as the Ctj used is noiser than that expected of a the full CMASS 
volume.). To guide the eye we plot the regions of constant a and e, as indicated in the legend (where the thicker line of each indicates 
the larger value). 



(e = 0, Q = 1), whereas if we would apply a geometric 
distoriton of |e|> 2% we should expect larger differences. 

In conclusion, the non-linear AP correction should be 
applied to avoid potential estimation biases. 



APPENDIX C: TESTING THE ALGORITHM 
ON HIGH S/N MOCKS 

We test our methodology by applying it on a set of 100 
mocks with higher S/N than those used in the final mock 
DR9 analysis. The motivation for this procedure is to sepa- 
rate between potential systematics and effects due to weak 
baryonic acoustic feature signals. 

The higher S/N mocks, called "stacked-mocks", are built 
by stacking the 600 PTHalo DR9- volume mocks by groups of 
six, providing us with one hundred realizations. For purposes 
of this analysis we divide the DR9 dj (see tj5.1.1|l by a factor 
of six. 

Figure [CTI shows distributions of (a||-{a||))/(Ta|| and 
(a_L-(ax))/f"a^ for the stacked mocks (top) and the DR9 
mocks (bottom) both pre- (left) and post-reconstruction 
(right). The quoted p- values are obtained when performing 
the standard Kolmogorov-Smirnov test between the distri- 
butions and a Gaussian one. 

We find that the stacked mock results yield various 
Gaussian (or symmetric) attributes not found in the DR9 
mock results. First, in the stacked mocks the means of the 
MCMC propositions are similar to the mode values, the 
standard deviations of the MCMC propositions are simi- 
lar to the 68CLr and they yield low skewness values of the 



marginalized ID likelihood distributions. As discussed in 
i]6.2l in the DR9- volume mocks we find large skewness caus- 
ing differences in these statistics. Using the DR9 mocks, we 
find in the that the modes and 68CLr are more reliable, as 
they are better defined. 

One of the most important Gaussian-like features found 
in the stacked-mock H'^'^jH, Dx/Dj^ results is that the 
scatter in the modes is similar to the mean of the uncer- 
tainties. This is not the case for the DR9-volume mocks, 
probably due to weak anisotropic baryonic acoustic feature 
detections. 

Finally, the stacked mock results (modes and uncertain- 
ties) are similar to those yielded when applying the same 
Cj^^ on the mock-mean signal (i.e., the mean signal of all 
600 mocks). We find this to be true for all eight combi- 
nations investigated: clustering wedges, multipoles; RPT- 
based, dewiggled templates; pre-, post-reconstruction. All 
results are presented in Table [CT] 

CI RPT-based vs. dewiggled templates 

As for preference of template (RPT-based vs. dewiggled) 
for constraining H^'^/H and Da/D^, when using the 
stacked mocks we find strong cross correlation coefficients 
of r ~ 0.9 — 1 in both modes and uncertainties. This com- 
parison shows no difference in uncertainties. The only odd- 
ity we find is that the dewiggled pre-reconstruction wedges 
and multipoles yield median (mean) biases of 1.4,0.9% 
(0.9, 1.0%) in H^'^/H modes, respectively, which is reduced 
post-reconstruction to 0.8,0.7% (0.7%). These H^'^/H bi- 
ases, when using the dewiggled model, do not appear when 



© 0000 RAS, MNRAS 000, 000-000 



BOSS DR9 H{z) and Da{z) from Clustering Wedges 23 
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Figure CI. (a||-(a||))/(Ta|| results (and similar for a±) of 100 6-stacked-mocks (Top) and 600 DR9 mocks (Bottom) pre- (Left) and 
post-reconstruction (Right). Results are for RPT-based clustering wedges. The p- values reflect K-S tests when comparing to a Gaussian 
distribution (blue lines). The p- values vary by template (RPT-based, dewiggled), and ^ statistic (clustering wedges, multipoles) used. 
The stacked- mocks yield p- values between 20 — 95%, where the DR9 mocks results have negligible p- values. 



Table CI. High S/N (6-stackcd) mock results 



^ (# of realizations) 



Actj I /o| I 



ax 



RPT-based wedges pre-rec (100) 1.002 ± 0.033 0.031 ± 0.006 0.996 ± 0.013 0.017 ±0.002 

RPT-based multipoles pre-rec (100) 1.002 ±0.033 0.030 ± 0.004 0.994 ± 0.013 0.016 ± 0.001 

RPT-based wedges post-Rec (100) 1.005 ±0.018 0.021 ± 0.003 0.996 ± 0.010 0.012 ± 0.002 

RPT-based multipoles post-rec (100) 1.003 ±0.016 0.022 ± 0.002 0.997 ±0.010 0.012 ± 0.001 

dewiggled wedges pre-rec (100) 1.014 ±0.034 0.032 ± 0.006 1.003 ±0.014 0.017 ±0.002 

dewiggled multipoles pre-rec (100) 1.009 ±0.032 0.029 ± 0.004 1.003 ±0.013 0.016 ± 0.001 

dewiggled wedges post-rec (100) 1.008 ±0.019 0.020 ± 0.003 1.000 ±0.011 0.012 ± 0.002 

dewiggled multipoles post-rec (100) 1.007 ±0.014 0.017 ± 0.001 1.001 ± 0.009 0.010 ± 0.001 



* The Oil and ax columns show the median and rms of the modes. 

* The Aa||/a|| and Aax/ax columns show the median and rms of the fractional uncertainties. 
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applied to the DR9-mocks. In those mocks, we find that the 
pre-reconstruction dewiggled model yields a bias of ~ 1% on 
determining Da/D^. 

In all four RPT-based cases (wedges, multipoles; pre-, 
post-recons truction) the mean biases of H^'^/H Da/D%^ 
are < 0.5%. ISanchez" et al.l (|2008h thoroughly analyze differ- 
ences between RPT-based and dewiggled and report that, 
when using the latter, one should expect systematic shifts 
in a due to the lack of a fc-mode coupling term. In i]5.2.2l we 
demonstrate that the post-reconstruction mocks do not pre- 
fer a template with Amc = 0, and hence suggest templates 
require a mode coupling term. 

For all the reasons above our choice of preference is the 
RPT-based template. 

C2 Clustering wedges vs. multipoles 

The stacked mocks show no significant difference regard- 
ing the constraining power of ^|| x and ^0,2 on H^'^/H or 
Da/D^; post-reconstruction RPT-based yields sub 0.1% 
differences. The cross correlation between the uncertainties 
of H^'^/H are found to be r ~ 0.6, 0.7 (dewiggled, RPT- 
based), and 0.88, 0.83 for Da/D*^. The pre-reconstruction 
templates yield similar results. 

We then ask if multipoles and wedges yield similar mode 
results. The post-reconstruction stacked mocks indicate r ~ 
0.80 for H^'^/H and r ~ 0.85 for Da/D^ in both RPT- 
based and dewiggled templates. Pre-reconstruction results 
yield similar correlations. 

For a visual of the results of the 3-stacked mocks, please 
refer to the bottom plot of Figure 1121 which is described in 

MM 

C3 Improvement due to reconstruction 

According to the stacked mock ^|| x (and hence also ^0,2), we 
find the uncertainty of H^'^/H improves by 32% and that 
for Da/D%'^ by 30%. 

The stacked mocks show that the H^'^/H modes should 
have a moderate correlation of r ~ 0.5 — 0.55 and Da/D^ 
of 0.5 — 0.6. For a visual of results from the 3-stacked mocks, 
please refer to the top plot of Figure 1121 which is described 
in El 

Another value of interest is the cross-correlation be- 
tween H^'^/H and Da/D^. With the stacked mocks 
we find this correlation to be of order r ~ —0.55 pre- 
reconstruction and r ~ —0.35 post-reconstruction. Also we 
find no correlation between a and e modes, as expected. 



